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Abstract 

We introduce and analyze several variants of a system of differential equations which model 
the dynamics of social outbursts, such as riots. The systems involve the coupling of an explicit 
variable representing the intensity of rioting activity and an underlying (implicit) field of social 
tension. Our models include the effects of exogenous and endogenous factors as well as various 
propagation mechanisms. From numerical and mathematical analysis of these models we show 
that the assumptions made on how different locations influence one another and how the tension 
in the system disperses play a major role on the qualitative behavior of bursts of social unrest. 
Furthermore, we analyze here various properties of these systems, such as the existence of trav¬ 
eling wave solutions, and formulate some new open mathematical problems which arise from our 
work. 


1 Introduction 

This article proposes a framework for describing the internal dynamics of riots, focusing on self¬ 
reinforcement and spatial diffusion mechanisms. The purpose here is not to explore the economics, 
social or political origins of riots, even less to discuss the legitimacy of any given riot. The approach 
in our work bears partial similarities to a recent literature on the use of mathematics in the analysis 
of uncivil and criminal activities |SDP+08[ iRNTnl IMSR+ 11 j . showing for instance that statistical 
regularities in crime patterns, together with insights from criminology, allow one to anticipate the 
evolution of such patterns, at least on short time scales. However, the phenomena and the models 
are quite different. We now first review some of the most common traits in riots to motivate some 
of the key ideas in our model. 

Civil disobedience and riots have been and continue to be means for populations or segments 
of populations to express their discontent towards their government or to react to certain events 
or political decisions |Moo78| . While these episodes of bursts of social activity are relatively rare 
compared to other phenomena, such as residential burglaries which essentially never cease, their 
current significance remains. If need be, the recent outbursts of protests and civil disobedience 
that broke out in Ferguson, Missouri (US) after the fatal shooting of Michael Brown by a police 
officer on August 9, 2014 reminded us of the relevance of this phenomenon |CL14| . Brown, who 
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was African-American, was short by Darren Wilson, a white officer, without apparent probable 
cause and this sparked a national and international conversation about the inadequate and tense 
relationship between the law enforcement and the community |LLB14j . A second shooting in St. 
Louis city reignited the protests and increased the tension between the African American community 
and the police force in Ferguson |Sall4| . 

Whether the tension between the population and their ruling government or the police arises 
from political events or decisions {e.g. the riots following the assassination of Julius Cesar in 44 BC 
in Rome, or the New York draft riots in 1863), new or increased taxes (such as the Moscow salt 
riots in 1648) , food scarcity |WS94| . high unemployment |Fla05| . police brutality |LLB14[ [BBJ12| . 
or racial tension, these events of social unrest are normally believed to have been triggered by a 
specific individual event. However, one can think of this “triggering event” as the the straw that 
broke the camel’s back. 

Consider, for example, the beating of Rodney G. King on March 3, 1991 by a group of policemen 
in Los Angeles, California (US) that was caught on video-tape and was made available for the world 
to see. Interestingly enough, it was not the beating of King that sparked the riots, but rather the 
injustice believed to be committed when the police officers involved were exonerated. Indeed, on 
April 29 of the next year the four police officers involved in this incident were acquitted and the 
first incident of the 1992 Los Angeles riots was reported only two hours later |Del95[ [Mnc09j . 

Of particular interest to us, is the case of the 2005 riots in France. The triggering event was 
the incident involving three young men who jumped into a power substation while being pursued 
by the police in Clichy-sous-Bois, one of the poorest suburbs of Paris, on October 27, 2005. Two 
of these young men died and this was the spark for the riots that spread throughout the country 
and lasted over three weeks |SVCB07] . As a final example, we mention the case of Mark Duggan 
who was shot on the chest by the police in Northern London on August 4, 2011, this triggered a 
four-day riot that spread throughout London |BBJ12j . 

These are only a few examples, but similar episodes continue to be observed throughout the 
world. Needless to say that the concept of civil disobedience is not new and has lead to many 
revolutions |Are72j . As Henry David Thoreau put it in his essay Civil Disobedience the idea behind 
civil disobedience is the belief that a just person must stand for what is right, which is not necessarily 
what is lawful... “// a thousand men were not to pay their tax hills this year, that would not be a 
violent and bloody measure, as it would be to pay them, and enable the State to commit violence and 
shed innocent blood. This is, in fact, the definition of a peaceable revolution, if any such is possible.^^ 

While it is widely accepted that all episodes of civil disobedience can be traced back to a single 
event, it is unclear what events are going to generate a cascade of civil unrest. In fact, there are 
many incidents that are very similar in nature to the triggering events mentioned above which do 
not generate bursts of rioting activity. For example, many incidents of shootings by police officers 
do not generate riots. If a community is content (employed, fed, educated, etc.) such incidents 
might be viewed as unfortunate occurrences, but citizens might not feel compelled to take to the 
streets and protest. However, if the social tension is sufficiently high, the shooting that might soon 
be forgotten in a content community, will have a high probability of igniting riots. In most serious 
riots in France, the triggering event was the death of a young individual in a poor neighborhood in 
situations where the police, or some other official authority, was involved. 

From the above examples we see that for a riot outburst to occur the social system needs to be 
“ripe,” in the sense that the social tension in the system needs to be sufficiently high. The social 
tension is a function of the economy, the police relationships with the community, the education 
level, as well as individual events. For example the visit of President Sarkozy to Clichy-sous-Bois 
on June 2005 where he stated that “we will clean up the city with a karcher,” likely raised the 
tension in this poor community. Also, in the recent years, solidarities based on religious issues in 
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poor neighborhood have started to play a role in generating local riots in France (as in Trappes, a 
suburb of Paris, July 2013). 

From our perspective, episodes of civil unrest require three factors: exogenous events (these 
include the triggering event, events that increase the tension in the system, and even “pre-triggering” 
events), endogenous factors (self-reinforcement in the system), and a sufficiently high social tension 
(a “ripe” system that is ready to experience self-reinforcement). Endogenous factors are those that 
are internal to a system, such as would be word-of-mouth type effects that promote the propagation 
and organization of the civil unrest. For example, during the French riots in 2005 the police 
announced on a daily basis the number of rioting events that occurred and the whole country was 
aware of the spread and level of the rioting activity. As a matter of fact, some believed that these 
announcements, if anything, were fueling the continuation of the riots. However, the government 
stopped announcing these numbers three days prior to the total cease of rioting activity and at 
this point the level of rioting was already quite small (less than a total of one hundred events per 
day). On the other hand, exogenous factors are external factors to the system that also affect the 
behavior of a system but in a different way. 

This dichotomy has been observed in various systems, such as world-wide-web searches and 
the number of times YouTube videos are viewed. Recently, in |MSB^11 Mohler and collaborators 


introduced the idea of modeling certain criminal activity, which experience repeat and near-repeat 
victimization [SDBT09] . as an epidemic-like phenomena using Hawke’s processes. Crane and Sor- 
nette in |CS08| introduced a method to determine the quality of highly viewed YouTube videos 
by using a Hawke’s process-like model to extract the effects of the exogenous factors versus the 
endogenous factors. As noted in |CS08| a great example of endogenous factors playing a significant 
role is the number of views of the Harry Potter trailer; on the other hand, the “tsunami” keyword 
search outburst was completely generated by the exogenous factor, mainly the tsunami that shook 
Japan in 2011. One of the advantages of the model introduced in |CS08j is that it affords the ability 
to extract the quality of videos. For example, a video with a sharp increase in the number of views 
and sudden decrease is more likely to have been boosted by an exogenous factor, but the quality 
of the video is probably not sufficiently good for viewers to pass along to others in their social 
network. On the other hand, the number of views of a video that experiences a slower but steadier 
increase is more likely to have been fueled by the quality of the video. An analogy can be drawn to 
civil unrest: episodes with a sudden spike and rapid self-relaxation are likely the mark of a strong 
exogenous factor and episodes which experience a slower but steady increase are likely to be fueled 
by the endogenous factors. The latter occurs when the system is experiencing a high and slowly 
decaying tension, and one can argue that these are more serious episodes or, at least, it is natural 
to expect that they will last longer. 

When dealing with systems that are inherently spatial a fourth factor is the inffuence that one 
location has on another. This effect has become extremely important due to the globalization of 
information brought about by the spread of technology and social media, a particularly suggestive 
example of the importance of this effect was observed in the “Arab Spring”, a revolutionary wave 
of riots, demonstrations, and protests that began in 2010 and spread throughout many countries 
including Tunisia (where the triggering event occurred), Egypt, Syria, and Libya - see |Lyn05| for 
a historical account and [T7m4] for a mathematical model related to revolutions. Therefore, it is 
of much interest to understand how the rioting activities spread spatially. What leads some riots 
to spread while others remain localized? In France, for example, there has been riots before and 
after 2005, but these riots remained at a local level - as in Vaulx-en-Velin, suburb of Lyon, October 
1992, the Sapins neighborhood in Rouen in January 1994, the La Duchere neighborhood in Lyon, 
October 1995, etc. In October 1995, in Vaulx-en-Velin again, a rather severe riot of about two 
hundred young people extended to the rest of the suburb of Lyon, but did not lead to riots on a 
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national scale. Similarly, in 2007 at Villiers-le-Bel, a two day riot propagated only to neighboring 
cities. A sociologist called these riots “les emeutes de la mort” [PerlOj . 

Early models of riots formation have focused on the emergence of a collective phenomena due 
to herding behaviors |Sch73[[Gra78j . Such models, formally related to Ising models in physics, have 
led to an important literature with various applications in social and economic sciences - for recent 
works, see [GNPSOQ], IBonl3j and references therein. In the recent years, the generic dynamics of 
riots and other social phenomena, as evoked above, has attracted more attention. In particular, 
several works have developed models based on the assumption that these systems are driven by 
Hawke’s processes |Haw - see for example |Oga98[ ICS08[ iLDD^ld) and references within. 

In this work we propose a model for the bursts of civil disobedience that includes the four 
mechanisms discussed above: endogenous factors, exogenous factors, sufficiently high social tension 
and influences, both local and non-local. In the bigger scope one of the ultimate objectives is to 
extract the strength of the exogenous factors versus the endogenous factors. 

As a first step, in this work we introduce a stochastic system on a network, which seems fitting for 
this application. Numerical realizations of the model illustrate a rich set of behavior of the system. 
We observe, especially in some parameter regimes, that the model behaves in a qualitatively similar 
way to what is observed in many real-world riots - capturing the global behavior without capturing 
the details. Moreover, we explore the effects that the non-local spread of information has on the 
spread of civil unrest. Through the development of the model we observe that it is necessary that 
the social-tension in the system also spreads in order for the system to experience a large scale 
burst of activity. To explore this phenomena mathematically, we then derive a system of nonlinear 
partial differential equations. 

This work presents a general model about riots and was initially motivated by the 2005 French 
riots. However, we believe that the family of models introduced here can be applied or adapted 
to many systems that experience bursts of activity followed by a period of relaxation. The general 
spirit of the model is to combine an explicit function which is observable with an implicit field, here 
the social tension. The explicit function - here the level of rioting activity, however it is defined 
- corresponds to actions that can be measured. The latter implicit field can be thought of as a 
potential field. The plausible existence of an underlying social field in collective social phenomena 
is stressed in |BBJ14| . We believe that this approach is relevant for a number of other situations 
where the introduction of such an underlying implicit field is warranted. For example, conflicts 
of various nature, with their escalation parts present some of the aspects we have described here. 
At the same time, models in the same vein but with bistable non-linearity in the implicit field are 
relevant to describe the loss of confidence, be it among people, towards organizations, media or 
towards the state of the economy, is such an instance where a buildup of distrust can enable a 
seemingly minor action or event to precipitate a complete loss of confidence. 

We also note that the system we introduce here has some similarities with models in neuroscience. 
Specifically it is related to models of neural dynamics that take into account synaptic depression 
|TPM98| . We give more details about this link in the next section (see discussion following equation 
(§). It would be interesting to further explore these analogies. 

We wish to emphasize that our objective here is to introduce simple models whose solutions 
exhibit the ‘stylized facts’ observed in riots, which can vary significantly from riot to riot. For 
example, the 2011 London riots was on the rise for four days and essentially ceased abruptly on 
day five, whereas the 2005 French riots took a period of about twenty-five days with long periods of 
increased activity and self-relaxation. There have been models for riots introduced in the literature 
previously - see for example |Bral2|. ILD14^ IDFWBl^ . which we discuss in more detail below. 
However, to the authors’ knowledge this is the first PDE model developed for riots, although 
the use of PDE systems to model urban crime has recently become an active field of research 
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SDP^OSl IBNlOj . A continuous model affords us the ability to prove rigorous spreading and decay 
estimates using PDE methods. It would be worthwhile, in future work, to move from a model 
describing the global behavior to one that actually captures more details about the actual spread 
of a given rioting activity. 

Lastly, we discuss some previous works that are especially related to our research. Theories on 
the role of contagion or diffusion processes in collective actions have been proposed already a long 
time ago (see e. g. LeBon 1895 |Bon05j . and references in |MyeOO| ), but formal mathematical models 
are quite recent. In |Bral2| . Braha introduced a non-linear spatial dynamical model for the global 
spread of civil unrest that includes short-range connections (describing geographic locations) and 
long-range connections (describing the effects of social networks and the media). In this work the 
author concludes that external causes such as those mentioned in the introduction (racial tension, 
food scarcity, etc.) are not necessary for the sudden outburst of civil disobedience. This is in contrast 
to our observation that the dynamics of these external causes are essential to fully understand any 
rioting activity. Davies and collaborators studied a model in |DFWB1^ , similar in nature to that of 
Braha, but which included the effects of police deterrence. This model was particularly concerned 
with data from the 2011 riots that took place in London. The authors tune the parameters of 
the model and obtain a simulated bursts of rioting activity that is qualitatively similar to what 
occurred during the 2011 riots. Finally, we mention a compartmentalized model for the dynamics 
of revolutions |LD14) . The authors of that work use a simple differential equation model that 
includes the effects of police repression and censorship. Under the assumption that for a revolution 
to propagate there needs to be a sufficient amount of protesting taking place, the authors divide 
the parameters of the model into regions that would lead to either stable or unstable regimes. 

While the models we introduce here contain similar ideas to those used in the works mentioned 
above, our model has a wider scope. It contributes the dynamics of the social tension in the 
system and a continuous model that allows us to rigorously analyze certain key characteristics of 
the spread of riots. The use of continuous models, which are less conventional in these contexts, 
to describe social phenomena has been popularized in recent years. Of particular interest are the 
models of Short and collaborators in |SDP^08| and of Berestycki and Nadal in |BN10| to describe 
the propagation of crime. These works introduce a notion of an invisible scalar field that measures 
the probability that a criminal activity occurs, this field is refereed to as the “attractiveness field” 
in |SDP^08| and the “willingness to commit a crime” in |BN10| . In the context of riots or civil 
unrest, this is analogous to a measure of the social tension. This concept is also found in the model 
of |Bral2| . which includes a measure of the political, social, and economic stress. 

Outline of the paper. In section we introduce the model on a network and illustrate the results 
of some numerical realizations of the system on a single site in section]^ We analyze the dynamics 
in the simplest case of a single site in section In the following section we perform and illustrate 
some numerical experiments on a network. We derive the continuous system in the case of local and 
non-local spread of the social tension in section We discuss the propagation of rioting activity in 
section 0 We conclude with a discussion in section 0 


2 Description of the model 

It is natural to consider a network of N nodes, where each node represents a location that is prone to 
rioting activity. These nodes correspond to the “urban clusters” of |Bral2j and they can represent, 
for example, cities in a country or neighborhoods within a city that are likely locations for the 
gatherings of people who are protesting or participating in more violent and destructive activities, 
such as arson or looting. Let us denote this network by A/”, we discuss the connections between 
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the nodes in a network shortly. We assume that for any node s ^ J\f there is an explicit field that 
measures the level of rioting aetivity. The level of activity is dynamic in time and we denote it by 
A(s,t). We further assume that there is a base intensity rate \h{s) which can vary between nodes 
on the network. This base level determines the low recurrent activity that occurs in the absence 
of any unusual factors. For example, according to media sources, the typical number of burnt cars 
in France is between fifty to one-hundred per night, with as many as three to four hundred burnt 
during special times, such as New Year’s eve. Interestingly, since around 1999 the burning of cars, 
in the absence of riots, has become a national sport in France. In fact, it is estimated that 10-30% 
are insurance crimes. This is clearly a component of the self-reinforcing mechanism: some people 
might burn their own cars after other cars have been burnt in their neighborhood so as to make the 
insurance company believe that they are victims of these events. 

Under the assumption that the number of rioting activities follows an inhomogeneous Poisson 
process, the X{s^t) would correspond to the intensity of the process, which is the expected number 
of rioting events. In this case, the expected number of events that occur during the time interval 
(a, h) at node s is given by: 

\{s^{a^h))— / X{s^t) dt. 

J a 

Assuming that the system is “ripe,” then the occurrence of a triggering event will spark a 
movement that gains momentum and is self-reinforcing. The endogenous effect (or self-exeitement) 
is built in the dynamics of A(s,t). Of course, there is a natural saturation limit, e.g. there is a 
maximal number of building and cars that can be destroyed. Thus, we begin with the following 
dynamics of the level of rioting activity: 

^A(s,t) = -w(A(s,t) - A6 (s)) + G{X{s,t)), (1) 

where cj is the natural mean reverting parameter of the level of rioting activity if there is no self¬ 
reinforcing activity. For simplicity in most of this paper we will assume A^ = 0. However, including 
a non zero base level of activity will be essential for future validation of the model with data. The 
self-reinforcement mechanism is modeled by the function G{z) which satisfies: 

G{z) > 0 for z G (0, zo), G(0) = 0 and G{z) < 0 ioi z > zq. (2) 

An example of this is a KPP-type term: G{z) = z{zo — z) for z G (0, zq) for some zq > 0. 

Since riots are bursts of social activity triggered by exogenous events, but not all external events 
(similar in nature) lead to riots, we assume that the systems must be “sufficiently ripe.” To express 
this mathematically, we introduce an implicit variable that represents the readiness of the system 
to experience these bursts. We refer to this scalar field as the soeial tension and denote it by 
a{s^t). Naturally, this value can vary from cluster to cluster and it is dynamic in time. Moreover, 
we assume that external events, such as, controversial political remarks, lack of social justice, high 
unemployment rates, police brutality, and the triggering events, tend to increase the tension in the 
system. In a way the social tension measures the level of resentment that a community or population 
feels toward the authority they are facing. Thus, we think of the kindling of riots as being analogous 
to flame propagation: the endogenous factors take substantial effect when the tension has reached 
a eritieal tension. This leads us to update 0 as follows: 

^A(s,t) = -Lo{X{s,t) - Xbis)) + r{a{s,t))G{X{s,t)), (3) 

where r{z) is, for example, a sigmoid function: 

~ 1 _|_ ^-I3(z-a) ’ 
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where /3 > 0 provides a measure of the transition slope between a relaxed state (non-excited state) 
and an exeited state. In other words, it provides a measure of how fast the transition is between 
a system that does not include the endogenous factors and a system with the full-force of these 
factors. The critical tension is denoted by a. Note, that in the limit as /3 approaches infinity 
r{a) approaches a step function: then as soon as the tension is above the critical threshold the 
endogenous factors are in full-force. We refer to Figure for an illustration of this transition 
function. The idea of a critical threshold was introduced in |LD14| . where the authors make the 
assumption that the number of protests has to be sufficiently large before it begins to grow into a 
revolution. This assumptions leads to a bistable ordinary differential equation. However, we note 
that the critical value is assumed directly in the level of rioting activity, which stands in contrast 
with our model. 

It is clear that the dynamics of a{s^t) in the process are crucial. Based on the above discussion 
we first assume that the exogenous factors increase the social tension in the system. It is noteworthy 
to mention that there are numerous factors in a society, which are present at all times and do not 
correspond to a single event, but that also significantly affect the level of tension that a community 
experiences. These factors include, but are not limited to, the state of the economy, unemployment 
rates, and political tensions. We include the exogenous factors that can be pin-pointed to a particu¬ 
lar time and place as “point sources” in our model. These occurrences are deterministic and should 
be clear indicators that the social tension will increase. In general, if n exogenous events occur 
at the times and the locations {si^tiY^ C A/* x {t > 0} , the source term produced by this effect is 
given by YY!i=i^i^t=ti,s=si^ where Ai measures the intensity of the exogenous events and 6t=ti,s=si 
is the Dirac delta centered at Si). To include the more constant factors, such as the state of the 
economy, we introduce a source term 0 ^ 5 ( 5 ). For example, if the economy is on a downward turn, 
the social tension will increase: 0 ^ 5 ( 5 ) > 0 during this period. The second important hypothesis 
is that, in absence of exogenous inputs, the social tension tends to decay (hence r{a) will decay). 
This hypothesis is somewhat analogous to Myers’ proposal that the riot “infectiousness” decays 
gradually over time |Mye00| , which he tests with an econometric approach in the case of the US 
racial riots in the 60s. Incorporating these effects into the model gives: 

d ^ 

—Q!(s, t) = Ai5t=ti,s=si - h{X)a{s, t) + 9ab{s). 

1=1 

The function h{\) represents the effect that riots or protests can have on the tension: the higher 
the level of rioting activity the slower the tension decays. For example, h{X) could have the form 


h{\) — 0 exp(—p A) or h{X) — 9 ( 1 -h 

Ai 


-p 


( 4 ) 


The parameter 9 = h{0) measures the natural decrease of the tension per unit of time, which sets 
the natural timescale over which the exogenous factors have an effect. The equation is set in such a 
way that in the absence of shocks (Ai = 0) and when there is no rioting (i.e. A = 0), then, the social 
tension field a reverts to the base rate a^. Indeed, the equation then reduces to d = —9(a — as). 

The parameters p and Ai control the influence that A has on the decay of a. In the following, 
we will set Ai = 1 and assume a slow decay, taking h(X) = 9/(1 -\- A)^ with 0 < p < 1. See Figure 
lb for an illustration of this function with two different values of p > 0. However, from the point of 
view of modeling, both the cases p > 0 and p < 0 make sense. They describe different situations. 
We will discuss the case of a fast decay of activity, that is p < 0, in further work. 

The timescale over which the exogenous factors have an effect can be different from the timescale 
over which the endogenous factors have an effect. For example, it took one year for the Los Angeles 
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riots to begin after the initial release of the video showing the beating of King. However, the tension 
had been building up and the exoneration of the policemen responsible for the beating increased 
the tension above the critical tension. On the other hand, the riots only lasted six days. However, 
these six days were very intense and the total number of deaths surpassed that of any other riot in 
the United States with the exception of the New York city draft riots of 1863 |Fiz06| . This leads 
us to the condition that oj > 9^ so that the exogenous effects can be observed over a longer period 
of time than the endogenous effects. 




(a) Transition function: r(z) (b) Decay control function: h{z) 

Figure 1: (a) Illustration of the function r{z) with /3 = 1 and /3 = 10. (b) Illustration of the function 
h{z) with 0 = 1 for p = 1 and p = 0.5. 

Including all these elements, the model then takes the form of the following system: 

[- Xbis)) + r{a{s,t))G{X{s,t)), 

1 ^ Ai5t=ti,s=si - h{X)a{s, t) + Oabis). 

This system of equations has some similarities with the ones of the neural field equations in the 
presence of synaptic depressions, even though the non-linearities are not of the same form in the 
two cases. These equations, derived from the classical Wilson - Cowan equations |WC72| . have the 
form of a mean field model describing neural dynamics. Tsodyks, Pawelzik and Markram |TPM98| 
have introduced an extension of this model that takes into account synaptic depression, that is, 
the decrease of signal transmission at the synapse when there is a high level of activity. One can 
think of the rioting activity level here as the analogous of the neural activity of a population of 
neurons, and the social tension as the analogous of the quantity accounting for synaptic depression. 
The latter represents the amount of resources available at the synapses for signal transmission. In 
both cases, the second quantity modulates the reinforcement mechanism in the activity dynamics, 
and relaxes with a time scale depending on the activity level. However, in addition to the non- 
linearities, there are two main differences. First, in the neural context, there is no direct external 
inputs to the synaptic depression field (hence no shock terms as above). Second, when modeling 
the spatio-temporal neural dynamics, synaptic depression acts on the inffuence between different 
neurons or locations |TPM98[|BK10| . whereas in our case, we assume the social tension to be a field 
acting specifically on the feedback from the local activity on itself but not on the inffuence from 



















other locations (see below). In other words, here, the social tension field enters as a modulation of 
the self-excitatory local dynamics of rioting activity. 

In the form (§, the model does not include any effects that one node might have on another. 
However, recent popularization of social media has enabled a kind of globalization of information 
and it is now rare to find isolated communities. Indeed, what happens in one city is quickly known 
world-wide and therefore it is natural to assume that the endogenous and exogenous effects spread in 
potentially non-local ways. The role of the use of television networks, cell-phones, and social media, 
to name a few, is the subject of various studies |Alcl li [GBBHRMTT] . In particular, the works |Bral21 
IDFWBl^ ILD14j . which we discussed in the introduction, either include or note the importance of 
the spatial component and diffusion of information due not only to geographical proximity but also 
to social connections. Furthermore, in |BBJ12j the authors analyze the data from the 2011 London 
riots and conclude that there is evidence for the diffusion and clustering of rioting activities, which 
is a call for spatial models that explore such factors. Different neighborhoods, towns, cities, or even 
countries are influenced in spatially heterogeneous ways. For example, Marseilles, France is reputed 
to be a self-centered city that is not influenced by the rest of France, whereas rural communities 
might be highly influenced by urban centers. We assume that the social tension in a node spreads 
to nodes which are within their social network. This communications or social connections are 
encoded in the matrix C = {cij)nxn with 

f 1 if node i has a social or communication connection to node j, 

[0 otherwise. 

At the same time, it is also natural to account for geographic proximity effects. One is naturally 
concerned about what happens to neighbors more so than what happens globally. Indeed, it is clearly 
observed that the level of rioting activity spreads locally from a place to neighboring districts or 
cities. To account for this we assume that the level of rioting activity diffuses to geographically 
neighboring locations. To quantify the geographic proximity let us define the matrix V — (%)nxn 
such that: 


_ J 1 if node i is a neighbor of node j, 

[0 otherwise. 

We denote the degree of a node s, i.e. the number of geographic neighbors (or edges) that node s has, 
by dv{s) and the number of social connections by dc{s). Let 77 represent the total influence that the 
nodes connected to s have on node s. The actual effect that node s' has on node s is proportional to 
the total number of nodes connected to node s. Thus, when we include the influence of neighboring 
nodes, the dynamics of the level of rioting activity is defined by: 

^A(s, t) = t) + r(Q!(s, t))G{\{s, t)) - cj(A(s, t) - Xb{s)), 

which can be written using the graph Laplacian A^: 

t) = j^AgX{s, t) + kX{s, t) + r{a{s, t))G{X{s, t)), (6) 

where k. — rj — uo. The graph here is associated to neighboring locations. The choice of the graph 
Laplacian is the simplest and most convenient for an initial analysis. However, it must be said that 
if A(x,t) represents, for example, the number of people (or fraction of the population) which are 
protesting or rioting, then more complex topologies must be included in order to move toward a 
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more realistic model. In particular, the inclusion of transportation networks will be crucial. These 
heterogeneities could lead to non-local, fractional, or non-linear diffusion. Non-linear diffusion 
arises naturally if one takes into account that the diffusion should really be proportional to the 
level of rioting activity. Of course, the behavior of a model with the more general diffusion can be 
significantly different; however, the analysis of these cases is beyond the scope of this paper. 

From equation Q we see that the natural decay of the level of criminal activity must be larger 
than the total inffuence of the neighbors of node s. This gives rise to condition that cj > ij to 
guarantee the eventual decay or rioting activity. We assume that the inffuence on the social tension 
is governed by the social/communications network (which could be non-local geographically). In¬ 
corporating this into the model gives the dynamics of the social tension is governed by the following 
equation: 


d ^ 

—Q!(s,f) = ——'^Csja{j,t) + '^Ai5t=o^s=si - h{X{s,t))a{s,t) + eab{s). (7) 

^ ^ j i=l 

Combining equations Q and Q yields the final system on the network. This system can easily be 
generalized to include weights on the inffuence between any two nodes, in which case, we would use 
the weighted graph Laplacian. The following table summarizes the parameters. 


Parameters 

Description 

(jj 

Decay rate of the rioting activity level A. 

h 

Base rioting activity level. 

Ai 

Strength of the shock at time ti and location Si. 

9 

Decay rate of the social tension value a. 

P 

Level of inffuence that A has on the decay of the social tension. 

Ub 

Base social tension value. 

/? 

Sharpness of the transition between the relaxed state and excited state. 

a 

Critical social tension value. 

V 

Strength of the inffuence of neighboring nodes. 


We can also consider a stochastic version of this model, which takes into account, for example, 
the effects of the media and the climate. Let Xf represent a Brownian or Levy process, then we 
obtain the stochastic version of the model: 


dX{s^ t) — - —A^A(5, t)dt -h /^A(5, t)dt + r{a{s^ t))G{X{s^ t))dt + aX{s^ t)dXt^ 
dv[s) 

n 

da{s, t) = -jXr Csja{j, t)dt + W ^i5t=o,s=si - (^(A(s, t))a{s, t) - 6ab{s)) dt. 
dc{s) ^ ^ 

7 1=1 


(8a) 

(8b) 


System Q is a coupled system of stochastic differential equations. 


3 Numerical experiments of the single site model 

To demonstrate the ffexibility of the model to capture the various global behaviors observed in a 
variety of real-world riots, we begin with some numerical experiments. For this purpose, we consider 
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the following model on single node without any noise and = 0: 


f\{t) = + r{a{t))G{X{t)), 

(9a) 

= f^Ai5t=ti - h{\{t))a{t). 

(9b) 

A(0) = Ao and (a(0) = cto, 

(9c) 


with either one or two shocks (n = 1 or n = 2). 

1. Slow relaxation of rioting activity: The first simulation results in an outburst of rioting activity 
that relaxes slowly - refer to Figure The outburst was a result of one shock at time t = 0 
that was intense enough to push the social tension above the critical threshold. This type of 
relaxation was observed during the 2005 French riots. 

2. Fast relaxation of the rioting activity: The second simulation results in an outburst of rioting 
activity that suddenly decreases - refer to Figure This type of relaxation was observed 
during the 2011 London riots. The driving parameter here is [3 (all other parameters where 
unchanged from the previous simulation), as the transition in the function r is sharper the 
self-relaxation (and self-excitation) is also sharper. 

3. Delayed outburst of activity: The third simulation illustrates the case of two exogenous events: 
the first event occurring at time t = 0 and the second event occurring at t = 12. The first 
external event is not strong enough to lead to an outburst of activity, but it does increase the 
tension in the system. Thus, when the second external event occurs, this drives the social 
tension above the critical threshold leading to what we call a delayed burst of activity - refer 
to Figure This is similar to what happened in the 2001 L.A. riots: the first event would 
correspond to the beating of King and the second event to the exoneration of the police officers 
involved. 


Two bursts of activity: The final simulation results in two bursts of rioting activity. In this 
case the first shock was strong enough to lead to a bursts of activity that settled down before 
a second shock occurred (of smaller intensity), which reignited the activity - see Figure 2d 
This was observed, for example, in the recent protests in Ferguson, Missouri. 


4 Analysis of the single site model 

We now provide a more rigorous analysis of the system given by Q by first analyzing it in the 
absence of shocks, and then considering the case of a single shock {n — 1), and finally the case of 
repeated shocks. 

4.1 Absence of shocks 

We begin with the study of the system: 


f\{t) = $(A(t),a(t)), $(A, a) 

:— —Lv{\ — A^) -h r(ct)G(A), 

(lOa) 

II 

(P 

:= 9as — ah{\) 

(10b) 

A(0) = Ao and ct(0) = cto, 


(10c) 


where is assumed to be small relative to the maximum of the function G. 
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Case with = 0 and = 0. A preliminary remark is that the dynamics dX/dt = $(A(t),Q;) 
for a fixed a value, has as stable fixed point at A = 0 for cj > r((a)G'(0). On the other hand, for 
(jj < r{a)G'{0)^ the stable fixed point is A*(q;) > 0 which is the non-zero solution of $(A,d) = 0 - 
for example A* = zq — uj/r{a) when G{z) = — z). We thus choose the parameters such that 

G'(0)r(0) < (jj < G'{0) and lim r{a) = 1, (11) 


so that, at constant cr, the attractive state is a no-riot state for small tension and a non-zero rioting 
state for large tension. 

The dynamics are easily understood by looking at the nullclines of (10) in the plane (cr. A). The 
nullcline associated to A, $(A, a) = 0, is given by A = 0, and, for a > ac with r{ac) — Vc uj/G\G)^ 
the function Ai((a), which increases from 0 at a = o^c to its maximum value A* < zq as a ^ oc 
{r{a) ^ 1 in this limit). This is illustrated on Figure]^ for G = — z). During the evolution of 

A increases everywhere under this curve (a, Ai) if A > 0, and decreases elsewhere. 

Since = 0, the nullcline associated to a is simply cr = 0: indeed da{t)/dt is strictly negative 
for any a > 0. One can then conclude that there is a single attractive fixed point, a = 0,A = 0. 
If the initial values are in the domain left/above the A-nullcline, this fixed point is reached with 
continuously decreasing a and A. If one starts below the nullcline, with X{t = 0) > 0, the fixed 
point is reached after an excursion at high A. While a is permanently decreasing, A increases until 
the nullcline is reached with a null slope dXjda = 0, and then decreases towards 0. If A = 0 and 
a > ttc, the system is at an unstable fixed point. However, if one adds any small perturbation, 
the system will be driven into the A > 0 domain where the rioting excursion will occur. For an 
illustration, see Figure]^ (the A—nullcline is the dashed-red curve). 


Case with > 0, A^ > 0. When A^ > 0, the fixed point at A = 0 is replaced by a fixed point 
at a small A value. Along the A—nullcline, as a increases, A remains almost constant up to o^c, and 
then A increases sharply as for the case of A^ = 0. Because A > 0, there is no more an unstable 
fixed point at a small A value. 

For 0^5 > 0, the a—nullcline is a curve close to the vertical line for small and moves towards 
the right with smaller slopes. For small the two nullclines intersect in the range where A is small, 
and this gives the unique stable fixed point: starting from any initial condition, the trajectory will 
end at this point, with a riot-excursion if the initial values (ao, Aq are on the right of the A—nulcline. 
For a large value, the intersection is in the range where A is large. The unique fixed point is 
with a high level of rioting activity. If (3 is larger than some threshold, there is an intermediate 
domain of a values for which the two nullclines have three intersects (see Figure [4b|). In that case 
the intermediate intersect is an unstable fixed point, and the two others are stable fixed points: 
one with low rioting activity, one with high rioting activity. One has here a structure with a 
classical discontinuous, order’ transition, leading to abrupt changes of behavior and hysteresis 
phenomena. If increases starting from a very small value, there is a critical value al at which the 
fixed point with high A value appears. As increases, the system may remains on the low fixed 
point, until a second critical value, where the lower fixed point disappears: the system suddenly 
jumps to the state of high rioting activity. If now decreases smoothly, the system will remain on 
the high A fixed until = al where the high fixed point disappears, and the system is abruptly 
brought back to the no riot state. This is reminiscent of the ending of the 2011 London riots, which 
stopped unexpectedly abruptly. The decrease in temperature and arrival of rain has been among 
the factor contributing to this event. One may assume that the weather contributes to the level of 
at. 
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(a) Slow self-relaxation 


(b) Sharp self-relaxation 




(c) Delayed burst of activity 


(d) Double burst of activity 


Figure 2: Simulations of the system given by|^. The solutions illustrated in 2a use parameters: 
zq = 10; cj = .2; a = 5; 0 = .1; p = 1;/3 = 10, a = 6,ti = 0,^2 = 12. The solutions illustrated in[^ 


use parameters: zq = 10;u; = .2;A = 6;d = .1; p = l;/3 = l,a = 6, ti = 0. The solutions illustrated 
misuse parameters: zq = 10; cj = .3; A = 8; = .3; p = 1, j8 = 100, a = 6.. Finally, the solutions 


illustrated in 2d use parameters: zq = 10; cj = .3; Ai = 6; A 2 = 3, 0 = A; p = 1,/3 = 1, a = 6. 
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Remark 1. From the Poincare-Bendixon theorem it follows that there is always convergence of 
the solutions to [10) to a stationary state. Indeed, the Poincare-Bendixon theorem states that a 
trajectory necessarily either converges to a stationary point or converges to a cycle. But here, owing 
to the second equation in which the right-hand side has a sign, there are no cycles. Note that this 
theorem holds for both cases p < 0 and p > 0 it is only the sign that matters. 


4.2 Adding a single shock 

First consider the case = 0 and = 0. With an initial condition at small a and A, a single shock 
A6t=o on da sends the system at a higher value of a. If the nullcline is not crossed, both a and 
A then decrease smoothly towards the fixed point (0,0). If the nullcline is crossed, the trajectory 
makes a rioting excursion as described above. If the initial value of A is zero, a strong shock sends 
the system onto the unstable fixed point A = 0. However, if one adds any small perturbation, the 
system will be driven in the A > 0 domain where the riot excursion will occur. 

In the case when the base level of rioting activity is zero and = 0 we observe an eventual 
self-relaxation. The formal proof of eventual relaxation is as follows. 

Proposition 1 (Burst of activity with eventual self-relaxation). Let n = 1 and any Ai > 0. If \{t) 
and a{t) are solutions to Q with G'(0)r(0) < u, and X{t = 0) > 0, then 

lim X(t) = 0 and lim a(t) = 0. 

t^oo t^oo 

Proof. We can solve for the social tension explicitly: 

a{t) = Aie-foh(Hs))ds^ 


which we then substitute into the equation for A to obtain: 

d\{t) = -ujX{t)dt + r ds)\ G{X{t))dt. 


Recalling that A* is the maximum value of A, we have the following lower bound 

I * - rh' 


which allows us to conclude. 


□ 


Even though the level of rioting activity eventually ceases (or relaxes to its base level of activity 
Ao) as does the social tension in the system, it is clear that the higher the intensity of the triggering 
event the longer the bursts of rioting activity will last. Indeed, if the shock is strong enough, the 
maximum value of A along its trajectory is reached on the nullcline in the domain where the nullcline 
is in its asymptotic regime A ^ A*, so that the activity remains at its maximum for a long period 
time before being able to decrease - see again Figure for an illustration. In the same vein, the 
parameter 6 plays, in some sense, an even more important role on how fast the solution decays. 
Indeed, a decays exponentially fast at a rate which depends on 0: smaller values of 9 lead to a 
slower decay. More formally, we prove that when the shock is sufficiently large the intensity will 
remain close to its maximum value for long periods of time. 

Proposition 2 (Long-periods of rioting activity due to a strong shock). Given arbitrary L > 0 and 
5 > 0, there exists Aq = Ao{6, L) and to > 0 such that if A > Aq then 

X{t) > A* — 5 V t G [to5 to ~\~ L\. 
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Proof. From the a equation (10a) we get 


a{t) = Ae--^oHHs))ds^ 

where A can be interpreted as the initial shock. To take into account a pre-existing base rate we 
could also consider A = A' + 9ai) where A' is the initial shock and the result is the same. Since h is 
bounded from above, from this formula we infer that there exists some r > 0 such that a{t) > Ae~^^ 
for all t >0. Therefore, given 77 > 0, T > 0, there exists Aq such that for all A> Aq we get 

>1 — 77 , for all t G [0,T]. 


Then, consider the equation 


y = -uiy + {l-r])G{y), y(0) = Aq > 0. 

Clearly, \{t) > y{t) for all t G [0,T]. By our assumptions, this system admits a unique globally 
asymptotically stable equilibrium > 0. Hence, given 5 > 0, there exists to > 0 such that 
y{t) > Zr^ — 512^ for all t > to- Now, ^ A* as 77 ^ 0 and therefore we can choose 77 sufficiently 
small so that Zr^ > A* — 5/2. Gathering these properties, we see that given 5 > 0 and L > 0, we 
can choose accordingly 77 > 0, then to > 0, and then T = to + T for which we get Aq such that 
for all A > Ao, A(t) > A* — 5 on the interval of time [to, to + L]. This completes the proof of the 
proposition. □ 

Figure illustrates the dynamics for A^ > 0,0^^ > 0. The general features described in the 
previous case remain, with two important aspects. First, as discussed previously, for larger than 
some threshold there is a fixed point at high value of A: obviously a strong enough shock will make 
the system evolve towards this high activity state. Maybe more interesting is the behavior below 
the transition. Suppose the system is ‘at rest’, that is at its fixed point in the absence of shock. 
Then any shock will drive the system in the domain where X{t) increases. However, if the shock is 
small enough, the system remains into the domain where the A—nullcline is almost flat, the rioting 
activity remains very small. See Figures |5a| and 

4.3 Repeated shocks 

If multiple exogenous events occur, their intensity and proximity in time will play a role on whether 
the system remains excited or eventually relaxes. With strong enough shocks at a high frequency, the 
shocks will maintain the system within the domain below the A-nullcline, and since in that regime 
A always increases, the rioting activity will remain close to its maximum value A*. With shocks at 
moderate frequency, a given shock occurs while the system is in the regime where A is decreasing, 
and the shock sends back the system under the nullcline. One has then a (regular or irregular) cyclic 
activity. With weak enough stimulations, in amplitude and/or frequency, the system will eventually 
relax. For this single site case, whatever the shocks may be. A* is the maximum maximorum possible 
value of A, whenever the initial condition is smaller than A* (which is of course the case with the 
hypothesis of no riots at t = 0). 

To understand the roles that both the time between events and the intensities of the shocks 
play in either leading an eventual relaxation to the base activity value or the perpetual self-excited 
system, we explore the case when there are periodic exogenous events, each with intensity A and 
period T. It is useful to define the point source term: 

00 

S{t) = Aj2^t=iT- (12) 

i=0 


5b 
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(c) Dynamics 


Figure 3: Simulations of the system given by Q with uo = — 2^ ^ — 1^9 — .7^p — .7 

and varying intensity of the triggering event: a small shock A = 0.4, a moderate one A = 4 and 
a strong one, A = 30. Figure 3a shows the nullcline dXjda = 0, and the trajectories for the small 
(no riot generated) and the moderate (generation of riots) shocks. Figure [3b| shows the case of the 
strong shock (note the difference in the range of a values). In all cases the initial values are chosen 
with small a and A values. The dashed lines indicate the jump in a value due to the shock. Figure 


3c shows the time dynamics of the rioting activity in the three cases. 
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(a) Small P 


(b) Large /3 


Figure 4: Nullclines of the system given by (10) with G{z) = z{zo — z) and cj = .4, 2:9 = 2, a = 


1,0 = .7^p = .7, A 5 = 0.05. Left, [4al /3 = 3, and right, [4b| (3 = 6. For each case, three cr-nullclines 
are shown: 9ai) = 0.2,0.41 and 0.6. For /3 = 6 , a discontinuous transition occurs. 


Thus, we consider system ^ with the point source term given by S{t) as defined in ( 12 ). 


Proposition 3 (Periodic excitable systems). Let X{t) be the solution to system ([^ with point souree 
term given by Then, for any given T and for any given 5 > 0, there exists an A* = 5, Aq) 

sueh that for all A > A* we have 

liminf X(t) > A* — 5. 

t^oo 


This is a straightforward consequence of Proposition]^ Figuregives a numerical illustration 
in the case of a periodic excitation as above, starting from a positive value of Aq. If the stimulations 
are too weak or too spread out (below some thresholds Ac{T) or l/Tc{A)), the activity relaxes to 
null activity with damped oscillations. Above the threshold, the system converges to a limit cycle. 
The proposition states that these limit cycles are in regions closer and closer to A = A* as A gets 
larger. As shown on Figure]^ this is the same at fixed A when T becomes small (high frequency). 

We conjecture a similar result in the case of a stochastic source term: 


i=0 


(13) 


where each {Ai,ti} are specified by a compound Poisson process. More precisely, the time between 
events, U^i — U, are specified by {N{t) : t > 0} , a Poisson process with a prescribed intensity 
Ap. Moreover, {Ai : i > 1} are independent and identically distributed random variables with a 
prescribed distribution. Finally, we assume that {Ai : i > 1} are independent of {N{t) : t > 0} . 

Open questions. As t ^ oc, we conjecture that there is a limit distribution for X{t). Now, 
assume that A^ = 0. Then, given Ap , when the E(A^) are bounded and E(A^) ^ 0, we believe 
that this distribution converges in probability to 0 , and that, when E(A^) ^ 00 , this distribution 
converges weakly (or in probability) to the constant A*. It seems to be an interesting question to 
understand whether there is actually a phase transition, depending on E(A^), from a long time limit 
of Xt going to zero in probability to a long time limit of Xt approaching a non-zero distribution. This 
is what we observe in the periodic deterministic case form numerical simulations. In the random 
case, we note that whatever the value E(A^) is, we always have 


liminf A^ = 0 , 

t^oo 


lim sup Xt = A* 

t^OG 
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(b) Moderate at, small shock, zoom 


X 



(c) Near the transition 

Figure 5: Dynamics in the phase plane of the system given by with G{z) = z{zo — z) and 
cj = .4, zo = 2,/3 = 6, a = 1, 0 = .7,p = .7, = 0.05. Top and 9ai) = 0.2. Starting from the 

fixed point of the system without shocks, three trajectories are shown after stiulating the system 
with shocks of different amplitudes: A = 0.4,2 and 4. For a small enough shock (5b, A = 0.4), 
the system remains where the A—nullcline is almost flat, the rioting activity remains very small. 
Bottom,!^ 9ai) = 0.41. In this case, there is two stable fixed points. Under a strong enough shock 
- here two examples, A = 2 and A = 4-, the systems, initially at the low fixed point, ends at the 
fixed point with a high A value. 
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(b) Medium frequency, phase plane 




(c) High frequency (d) High frequency, zoom 

Figure 6: Simulations of the system given by Q with cj = .4, zq = 2,/3 = 3, a = 1, 0 = 0.7^p = .7 
under periodic stimulations at different mean frequencies with common amplitude A = 2. 


with probability one. Similar questions arise when E(A^) > 0 is given and one examines the 
inffuence of the intensity Ap, which represents here a frequency. We can even assume that the A^’s 
are constant. At high frequency, we expect the system to be close to the maximum value most 
of the time, while for low frequency, we expect it to be close to 0. Whether there is or for what 
distributions there is a phase transition is an open question. 

Some aspects of these questions are supported by numerical simulations. The dynamics with 
repeated shocks is illustrated on Figure for shocks of constant amplitude Ai = A but occurring 
at random times, with a mean frequency u — E(t^). As one would expect, there is a critical value 
z/c(A) such that for u < i^ciA) the systems eventually relax, and for z/ > i^dA) the system is in a 
sustained regime of bursts of excitations. It is observed that the activity is frequently close to the 
maximum A*, even at moderate rate. The fraction of time spent near A* increases as the mean 
frequency increases. 

5 Numerical experiments on a network 

To investigate the effects that spatial dispersal of information has on the spread of riots we illustrate 
and discuss some numerical realizations of the system defined in Q-Q in this section. We perform 
simulations on a network of one hundred urban centers that are on a square grid. The square grid 
gives the geographic neighbors of each node. 

1. Double-threshold phenomena: We study the case when one urban center is the only inffuence 
for all other centers. For example, we can think of this center as being Paris, which undoubt¬ 
edly has an inffuence on all of France. In order to explore the effect that the intensity of the 
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Figure 7: Simulations of the system given by Q with uo — — 2^ P — — 1^9 — 0.7,2? = 0.7 

under multi shocks stimulations at different mean frequencies with common amplitude A. The 
rioting activity, A(t)/A*, is shown as a function of time. 


triggering event has on the spread of the rioting activity, we perform a series of simulations 
with fixed the parameters uj — .2^9 — .3, zq = 10, /3 = 1, a = 100, p = .7, p = .2 and vary the 
strength of the triggering event A. In this case, we observe a double threshold phenomenon^ 
which we summarize below: 


(a) For A small (the simulation illustrates the case when A = 2) we observe that the rioting 
activity does not spread beyond the location of where the triggering event occurs. In fact, 
the intensity is not sufficiently high to provoke the spatial spread of rioting activity. Refer 
to Figures |8a||8^ for three different snapshots in time of this simulation. Red corresponds 
to high levels of activity and dark blue to zero levels of rioting activity. 

(b) For A of intermediate value (the simulation illustrates the case when A = 6) we observe 
that the rioting activity spreads in a local fashion: the burst of rioting activity affects 
nearest neighbors. It appears that the local diffusion of A(5,t) dominates in this regime, 
refer to Figures |8d}|8f| for the corresponding snapshots in time of this simulation. 

(c) For A large (the simulations illustrate the case when A = 10) the spread of rioting 
activity is non-local: the bursts of social activity jump infecting urban areas in a non¬ 
local manner. It appears that the non-local diffusion of a{sA) dominates in this regime, 
refer to Figures |8g||8l| 

The apparent threshold between the very localized rioting activity and the activity which 
spreads is not surprising, but the second observed threshold is. Indeed, from these numerical 
experiments we eonjeeture that there is a critical threshold value for the intensity, denote it 
by A*, such that if A < A* then the rioting activity spreads locally (the geographic network 
dominates in this regime) and if A > A* the rioting activity spreads non-locally (the social- 
network dominates in this regime). 


2. Delay effeet Consider a social network made up of two inffuential urban centers, for example, 
this could be Marseilles (M) and Paris (P). We assume that the triggering event occurs at 
time t = 0 in (P). If the initial triggering event is not sufficiently strong, we observe a nominal 
spread of the rioting activity and in fact the activity ceases relatively fast. Refer to Figure 


9b for an illustration. On the contrary, when a second exogenous event (even of significantly 


less strength) occurs at the second urban center (M) at a much later time {t = 30 for our 
simulation), we observe a much stronger and faster spread of activity both from (M) and (P). 
After the second exogenous event one observes a local spread of rioting activity from both 
(M) and (P). Refer to 9cpd for an illustration of this phenomena. Of course, this is due to 
the self-reinforcement between the two center (P) and (M), which leads to a delayed burst of 
activity. 
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6 Spatially continuous deterministic model 


While we can learn much from the discrete system and its numerical simulations, it is useful to 
work with the continuum limit in order to be able to perform more rigorous analysis. Although 
we recognize that a discrete model is more realistic for this application, the continuum models are 
useful in that they allow us to determine the intrinsic parameters of the system as well as to obtain 
estimates for the duration of an episode of rioting activity. Thus, we are motivated to introduce a 
couple of continuum versions of our models in this section. 

6.1 Local influence on the the social tension 

First, we consider the case when the diffusion of both the level of rioting activity and of the social 
tension is only due to geographic connections. This corresponds to the case when the matrix 
C — V. Let C (or as a first step in exploring this phenomena) represent our domain of 
interest. We are interested in the limit as the number of nodes approaches infinity {N oo): A(x, t) 
and a{x^t) are then defined for all x G li. We discretize and obtain a lattice and denote it by 
ri/ := “ Vi- nodes {xi.yj) e = Af 

correspond to the discrete locations where the rioting activity will take place, e.g. ‘urban clusters’ 
in the model introduced in |Bral2j . We are interested in the limit of the distance between the urban 
clusters approaching zero, i.e. Ax 0. For simplicity, we consider the case when the stochastic 
effects are negligible, which corresponds to the system given by equations (§ and 0. As the 
derivation of the continuum limit as both Ax and At go to zero is standard, we only mention that 
one uses the discrete Laplacian operator in order to keep track of mesh distance and take the limit 
in such a way that D ^ remains constant. Then, the corresponding limiting equation are 

the following: 


d 

dt 

d 

dt 


A(x, t) — DAX{x, t) + r((a(s, t))G(A(x, t)) — n\{x^ t), (14a) 

n 

a{x, t) = DAa{x, t) + Ai^ St=Q,s=si - (HMx, t)) - ri)a{x, t) + 9ab, (14b) 


1=1 


where hi — uo — rj. We always assume that /^ > 0. The existence of solutions to the Cauchy problem 
defined by (14) falls under classical theory. 


6.2 Single triggering event: Eventual decay of the rioting activity 

Let us consider the system with the following initial conditions: 

A(x,0) = Ao and (a(x,0) = 0. (15) 


We show that in the case when = 0 the bursts of social activity eventually cease though out the 
domain fl. That is, the we prove that the mass of the social tension and the level of rioting activity 
eventually approaches zero. 


Proposition 4 (Decay of mass). Let A(x,t) and a{x^t) he solutions to (14) with initial eonditions 
given by (15) for parameters sueh that 
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There exists fci, A ;2 > 0 sueh that 


||a(0, x) + 

1=1 

n 


( 16 ) 


i=l 


Furthermore, for any e > 0 there exists a > 0 sueh that A(t) < e for all t > T^. 
Proof. Let us first compute the dynamics of the mass of a, 


/ a{x,t) dx — Ai 2 _^ / ^t=n,x=xi dx — {h{\{x,t)) — rf) / a{x,t) dx. 

Jn Jn ’ Jn 


d 

dt jQ 


For simplicity let us denote y{t) = f^a(x,t) dx and k 2 — 0/[1 + X'Y — rj and ki = 9 — r]. Then we 
obtain the following upper and lower bounds: 

Th 71/ 

^ ^ / df=Q^x=xi dx kiy{t^ ^ ^ ^ / df=Q^x=xi dx k2y{tf 

Jn dt 

This gives the upper bound on the mass: 

1 rt ^ 

ds 

^ i=i 


<||a(0, x)||ie-'=2* + ^iY^ 


i=l 


Similarly, we obtain that 


y{t) > \\aiO,x)\\ie-'^^^ + 


i=l 


Now, let us consider the dynamics of the mass of A : 

^ f X{x,t) dx = f r{a)G{\) dx — n f \{x,t) dx 
dt Jq Jq Jq 

< r{g{t)) |fi| A* — / \{x,t) dx, 

Jn 

where, 

n 

g{t) = + 

i=l 

From this, we obtain the estimate 

ll^(a;>^)lli = ll^(a;,0)||ie“'"* + c [ r{g{s))e~'^^ ds. 

Jo 

Note that lima;^oo^(^) = 0 and this implies that lim^^oo'^(^(^)) = 0. Thus, for any e > 0 there 
exists a time T* = T*(t^, Ai,k 2 ^ a, /3, e) such that for all t > T"" A(t, x) < e. □ 
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6.3 Non-local influence on the tension 

Similarly, we can obtain a continuum model for the case when the social tension diffuses non-locally 
through the communications network C. In this case, we obtain a system with an integral operator. 


\t — DAX — {(jj — r])X + r(a)G(X) + cuX^ (17a) 

= r \ ^ / J{x,y)a{y,t) dy -h{X)a^eai,^ (17b) 

dy Jq 

The interaction potential J7(x, y) is either equal to zero or one, which is due to the choice of the 
social network (see the definition of the matrix C in section This, of course, can be generalized 
to include different weights. As an interpretation of J7(x, y) we can think that each location, x G 
has a domain of influence^ which is the set of locations which inffuence what happens in x. At the 
same time, each location inffuences some locations and not others, we refer to this as the range of 
influence of location x. Observe that the integral operator does not model diffusion and does not 
have to be symmetric. 

7 Traveling waves solutions 

One of the most important characteristics about the spread of rioting activity is the speed and the 
exact manner in which it spread. In fact, one can think of the spread of rioting activity as a front 
of high levels activity that is invading regions with the base level of activity. It is therefore natural 
to first look for the existence of traveling wave solutions. 

7.1 Local influence 

For simplicity, let us consider the model with local diffusion for both variables and with only one 
exogenous event. In fact, the shock can be included in the initial condition. For this purpose, 
consider the system: 


— A(x, t) — AA(x, t) + $((a(x, t), A(x, t)), 

^(a(x, t) — Aa{x^ t) + T((a(x, t), A(x, t)), 
A(x, 0) = Ao and a{x, 0) = Ai6x=x^ 


(18a) 

(18b) 

(18c) 


where $((a. A) := r{a)G{X) — nX and T((a, A) := —h{X)a -\- ya -\- 9ai). Solving da/dt = 0 gives the 
relationship: 


a{X) — 


6as 

h{X)-r]' 


(19) 
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Level of rioting activity at time: 4.375 


Level of rioting activity at time: 25.625 


Level of rioting activity at time: 50 






6 


i 


(a) A = 2, f ~ 4.3 


(h) A = 2, 25 


(c) A = 2, t = 50 


Level of rioting activity at time: 10.1563 




Level of rioting activity at time: 28.5938 



(d) ^ = 6, f ~ 10 


(e) A = 6, f ~ 15 


(t) A = 6, t = 50 



Figure 8: Simulations of the the network-dependent model (system given by (§-0) with (jj — 
.2,0 = . 3 , 2:0 = 10,a = 100,(p = .7, 77 = .2 and varying intensity of the triggering event. Figures 
8a||8c illustrate a simulation with a triggering event occurring at time zero with intensity A = 2. 
On the other hand. Figures |8d}|8f| illustrate a simulation where the triggering event has intensity 


^4 = 4. Figures SgHSl illustrate a simulation with an triggering event of intensity A = 10. 
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Level of rioting activity at time: 25 
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(a) t = 25 


(b) t = 50. 



(c) t — 19. 


(d) t — 45. 


Figure 9: Simulation with the social network with two influential nodes. Figures 9a]{% illustrate a 
single exogenous event with intensity A— Figures |9cl|M| illustrate the effect of a strong triggering 
event in (P) (red square in Fig[^ and a much waker second event A = 2 at t = 30, which leads to 
a delay spread in the rioting activity in (P). 


Substituting (19) into the equation $((a, A) = 0 then gives two steady state solutions: one 
corresponding to a non-excited state^ (<^i,Ai) := , and one corresponding to an excited 

state ((^ 2 , A 2 ) := A 2 ^ , where A 2 = A 2 (G, r). The excited state will always be stable for our 

choice of functions r{z) and G{z). However, if 


r((ai)G^(Ai) + 77 > /i(Ai) + fc. 


( 20 ) 


then the non-excited state will be unstable. In fact, as the critical tension decreases (with all 


other parameters fixed) the system defined by (18) goes from a regime where the two constant 
steady-states are stable into a regime where the non-excited steady-state becomes unstable. Refer 
to Figure [To| for an illustration of the reaction terms $((a(A), A) in the bistable regime (Figure [TOa] ) 
and monostable regime (Figure [TOb ). 

In fact, for a certain parameter regime there exist traveling wave solutions, ( 7 /^, 0, c) with c G M^, 
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such that for all z = x — ct the following holds: 


{ -ip^iz) + C'ip'{z) + $(</), 'Ip) = 0, 

C(p'{z) + C(p{z) + ^{(p,'ip) = 0, , . 

0<{i^{z),cP{z))<ia2,X2), ^ ^ 

V’(-oo) = 02, ^(-oo) = A2, ipi+oo) = Oi, (/>(+oo) = 0. 

The interesting observation is that the critical tension parameter a can determine if traveling wave 
solutions exist and whether they are unique. 




(a) Bistable system: a = 5 (b) Monostable system: a = 1 

Figure 10: Illustration of the $((a(A), A) with parameters: zq = 10,6c; = .2,0 = 0.05,77 = .01,p = .5. 


Theorem 1 (Critical threshold for traveling wave solutions). Let the he chosen so 

that ( [T^ has two steady-state solutions ((ai,0) and (0^2, A2) and such that (19) remains positive for 
all A G [0, A2]. There exists 0 < a^ < 00 so that: 


1. if a^ < a there exists a unique {f;{z)^(j){z)^c'' 
Furthermore, c* > 0. 


with 


up to translations, to {21), 


2. For 0 < a < a^ there exists a c* G and solutions , (j)c{z)) with z = x — ct to (21) for 

any c > c*. Furthermore, when c < c*, such waves do not exist. 

In either case, it holds that ip'{z) < 0 and (j)'{z) < 0 for all \z\ < 00 . 

The proof of Theorem is a direct application of two theorems from |VV94j . The only observa¬ 
tion needed to be made is that for fixed parameters u, 0, p,p, zq and k as hypothesized in Theorem 
([^, there exists a critical value a* such that if a > a* the (18) has two stable steady-state solutions. 


Thus, an application of Theorem 2.1 give part 1. On the other hand, for 0 < a < a* the non-excited 
state is unstable and an application of Theorem 2.2 gives 2. 

This result shows that certain initial data or triggering events can facilitate the propagation of 
rioting activity. Moreover, we observe a critical threshold develop for the critical social tension a. 
That is, for moderate levels of a the traveling wave can only move at a unique speed. However, if 
the a is sufficiently small then the traveling waves which facilitate the dispersal of rioting activity 
can spread at very large speeds. For example, a very intense triggering event could lead to a faster 
spread of the riots. We leave the study of the various qualitative properties of the traveling wave 
solutions from Theorem such as the asymptotic decay rates, for future work. 
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7.2 Numerical experiments on spreading 


The existence of traveling wave solutions can be useful in exploring the spreading speed of the riots 
and helps provide a global picture. However, the actual spread of the rioting activity in the 2005 
French data seems to have some very interesting properties, which would lead us to believe that 
the true solutions are not true traveling waves, although they do exhibit similar characteristics. To 
better explain this, let us denote by p{x) the peak number of events at location x and by t{x) the 
day in which the rioting activity peaked. In preliminary analysis of the spread of rioting activity 
during the 2005 French riots, we observe the following general trend: if the triggering event occurred 
in location x' at time t = 0 then the peak number of events decreased and was delayed in locations 
with the distance away from the rioting activity, in other words, p{y) is a decreasing function of 
\y — x'\ and t{y) is an increasing function of \y — x'\. This was only a general global trend and there 
were exceptions of course. In particular, large cities erupted first in some cases. We are able to 
reproduce this observation with the model 


Traveling wave-like solution: Fignres [T2a)| 1 2d| illustrate the numerical solutions where A(x, 0) = 
g-iOa: ^ triggering event occurring at location x = 0. The result is a wave-like solution 


whose crest decreases as it moves. In Figure [T^ one observes that the rioting activity grows 


and at time t = 1 the solution has a wave-like profile which connects the maximum level of 
rioting activity and the base level of rioting activity. At time t — 2 the solution continues 
to have a wave-like profile, however, it now connects a lower level of rioting activity to the 
zero level of rioting activity. Of course, this is not a true traveling wave as we see the crest 
decrease and the transition between the lowest level of activity and the highest level widen. 
We extract and illustrate from the evolution of this solution, the level of rioting activity at 
four different locations (x = 0,1, 2, 3) in Figure 11a Recall, that the triggering event occurred 
in location x = 0 and here the peak of the level of activity is the highest of the four and occurs 
the earliest. From there we observe that the peaks at each corresponding location decreases 
with the distance from the triggering event location. 


2. Spreading solution: Figures 13a||13d illustrate the numerical solutions where A(x,0) = Aq = 2 
and a triggering event occurring at location x = b. The result is a bump-like solution whose 
center is located at the triggering event location. Initially, this bump-like solution increases. 
In Figure 13a| one observed that the rioting activity grows and at time t = 1 it begins to 
simultaneously spread and decay. We extract and illustrate from the evolution of this solution, 
the level of rioting activity at four different locations {x — 2, 3,4, 5) in Figure lib Recall, 


that the triggering event occurred in location x — f) and here the peak of the level of activity 
is the highest of the four and occurs the earliest. From there we observe that the peaks at 
each corresponding location decreases with the distance from the triggering event location. 


7.3 Mixed local-non-local model 

The inclusion of social connections in the network leads to a (spatially) non-local spread of infor¬ 
mation in the social tension. To the authors’ knowledge this poses a new mathematical system for 
which there is much theory that needs to be developed. For example, while the global existence 
of solutions to the Cauchy problem defined by system (14) falls under classical theory, we are not 
aware of any proof of the global existence of solutions to (17) and hope to address this issue in the 
near future. Moreover, in the special case where there is a symmetric influence between nodes and 


27 




















this influence deceases with geographic distance we obtain a system of the following form: 


(22a) 

(22b) 


At = DAX — {lo — r])X + r{a)G{X) + uXb 

at ^ fj J{x - y)a{y,t) dy - - {h{X) + fj)a + 9ab + ^5o,s- 

The existence of traveling wave solutions in this system is an open problem. Although we conjecture 
that such solutions exists. 


8 Conclusion 


Our goal here is to introduce simple models for the dynamics of riots that include what we believe 
to be the essential ingredients to qualitatively capture the initiation, self-excitation, spreading, and 
relaxation of rioting activities. Without any of these factors, we believe that important qualitative 
information will be lost. In addition, we have introduced a way to model the social or global 
connections both on a network and on a continuum domain. The former is useful for the purpose 
of fitting data and the latter to extract some characteristics of the growth, spreading, and decay 
of the rioting activities. These models are mainly exploratory tools that can be used to gain some 
insight into the effect that the different mechanisms have on the rioting activity. 

One interesting outcome of this model is the double threshold phenomena that comes from the 
intensity of the triggering event. It was surprising to us that if the intensity of an event was 
sufficiently high the spreading of riots changed from being a local one to a non-local one. Of course, 
since this is an observation from numerical simulations it raises the question to verify this with 
rigorous analysis: it is an open problem to prove that such a threshold exists. 

While all rioting activities are similar in nature there are many contrasts between them as well 
- see for examole jNew m- We aim in this paper to construct a model that has the flexibility to 
describe various types of bursts of rioting activity. It is hoped that it could be used in the future 
to determine if, for example, a riot was exogenously driven or endogenously driven. This is another 
interesting direction of research that is worth pursuing. Likewise, the model allows for different 
possible regimes of influence of activity on the social tension. Here we have chosen to consider the 
effect when high activity slows down the relaxation of this tension. This is reflected in the choice 
p > 0 for the function h{X). In other situations, a high activity may actually generate a rapid 
decay of the social tension as some observations seem to suggest. This leads to choose p < 0 in this 
function. Also note that another possible choice for the evolution of the social tension field is 


-—a 

dt 


{s,t) = 


i=l 

We will discuss the case p < 0 and this 


Ai5t=ti,s=si - h{X){a{s,t) - ab{s)). 
variant elsewhere. 


(23) 
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(a) Level of rioting activity in time. 


(b) Level of rioting activity in time. 


Figure 11: On the left we observe the total number of rioting activities for four different locations 
with parameters and initial conditions corresponding to Figures fl2a|| 1 2d[ On the left we observe the 
total number of rioting activities for four different locations with parameters and initial conditions 
corresponding to Figures 13a[|13d 


Another interesting aspect of a riot that remains to be studied in detail is the process of “self¬ 
relaxation.” In fact, the reasons why riots end is recognized as quite unclear. For example, in the 
2005 French riots, riots were already strongly declining when the government adopted the strongest 
measures (notably an evening curfew). Rioters from the 2011 London riots reported to have stopped 
due to boredom, lack of targets, bad weather, fear of the police, and a call for peace by the father of 
a person killed during the riots |TLCllj . In the model, whenever there are no more shocks or other 
external inputs to the social tension, the rioting activity eventually ends (provided the background 
social tension is small enough). This is a key aspect of the model. The self-reinforcement effect 
leads to developing riots even though the social tension is already decreasing. However, when this 
social tension becomes small enough, it brings back the rioting activity to its baseline level. 

Some of the empirical factors mentioned in the introduction are clearly ingredients contributing 
to the decay of the social tension. It is well known that the weather plays a role in riots (everything 
else being equal, the likeliness of riots increases with the temperature). This factor can be thought 
of as modulating the baseline value of the social tension. Other empirical elements could be incor¬ 
porated into our model. Specifically, the effects of the police department or announcements for the 
police that the army will become involved (such as what happened in the 2011 London riots) can 
be incorporated as a negative shock. This is obviously the same for the call for peace - but why 
a specific call, or such an announcement by the police, are actually perceived as events of strong 
amplitude is a deep issue in sociology, outside the scope of this modeling approach. In this model, 
a negative shock on the social tension can bring the system into the regime of decreasing rioting 
activity. We would like to emphasize that a shock that would act directly on and only on the rioting 
activity would be inefficient - the system remaining in the regime where the activity increases. This 
bears some interesting consequences on what type of measures are more likely to ease a situation 
with severe unrests. A more detailed study of the effect of negative shocks is also a natural direction 
of future work. 

As already said, the purpose of this work is to establish models that characterize the global 
behavior of such systems. However, the ultimate objective is to use data from riots to validate 
this model. One perspective is the use of the 2011 London riots data which are available on line. 
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The case of the 2005 French riots would also be very interesting to study: the fact that the riots 
spread throughout the country makes this a particularly interesting data set to study the effects of 
non-local dispersal of the riots though the study of the proper social connections included in the 
model. In terms of stylized facts, our preliminary analysis seems to support that the model we 
propose here is coherent with the observations. 




(a) Solutions at t = 1. 


(b) Solutions at t = 2. 




(c) Solution at t = 3. 


(d) Solution at t = 4. 


Figure 12: Figures 12a||l2d illustrate the dynamics of the numerical solution to system (14) with 
exponential initial intensity at in the triggering event location x = 0 and with a triggering event 
intensity A = 50. The horizontal axis represents space. 


Finally, it is worth noting that modeling this social phenomena leads to the construction of 
new mathematical models, for which there is much mathematical theory to be developed. We hope 
that this article will encourage future work in this direction. In particular we emphasize the open 
problems related to the properties of the model with stochastic shocks and to the study of mixed local 
/non-local diffusion models. For the latter, solving the Cauchy problem, the construction of traveling 
fronts and the determination of the asymptotic speed of propagation are new types of problems which 
could lead to interesting mathematical developments. Lastly, the different geometries of a network 
and heterogeneities in a system will lead to varying types of diffusion. It would be an interesting 
problem to see how non-linear or fractional diffusion inffuence the the existence of traveling wave 
solutions and their propagation speed (if they exists). 
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(a) Solutions at t = .5. 


(b) Solutions at t — 1.5. 




(c) Solution at t = 2. 


(d) Solution at t = 2.5. 


Figure 13: Figures 13a 13d illustrate the dynamics of the numerical solution to system (14) with a 
constant initial conditions A(x,0) = 2 and a triggering event with intensity A = 100 occurring at 
location x = 0. The following parameters were used for both numerical experiments: zq = 10, a = 
100,6c; — .2^9 — 0.05,77 = .198. The horizontal axis represents space. 
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